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ABSTRACT 


A general review on the theory of plasma 
instabilities in a discharge tube is given for 
the helical type instabilities (azimuthal wave 
numbers m=1,2). The major effect causing the 
steady state radial density profile, as calcu- 
lated by previous authors, has been the boundary. 
They have considered the walls as an infinite sink, 
thereby causing a density gradient. After a sug- 
gestion by Glicksman we consider the self induced 
azimuthal magnetic field as the major cause of the 
radial density profile. The effect of the walls 
has been reduced by demanding that the tube radius 


be large. 
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CHAPTER I 


INTRODUCTION 


The motion of charged particles in the positive 
column of a discharge tube has been studied both the- 
oretically and experimentally in great detail by 
applying various configurations of electric and mag- 
netic fields. 

If a plasma is contained in a long cylindrical 
tube in which the walls are considered to be places 
of recombination, the electrons due to their lighter 
mass will tend to diffuse to the walls at a much 
faster rate than the ions. However, a charge separa- 
tion will result, and this charge separation will set 
up anelectric field which will retard the electron 
diffusion and at the same time enhance the diffusion 
of the positive ions. The net result is that the 
electrons and positive ions will tend to diffuse 
together at an intermediate rate. This is called 
ambipolar diffusion. If a constant electric field 
is applied parallel to the axis of the tube the diffu- 
Sion perpendicular to the electric field will remain 
ambipolar. However, the motion of the particles along 
the electric field will be controlled by the electric 
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The motion of a charged particle in a constant 
Magnetic field, with a constant electric field applied 
parallel to the magnetic field, may be described as 
the sum of a gyration around the magnetic field lines, 
a drift motion across the lines and a motion parallel 
to the field. In an ionized gas (plasma) this motion 
of particles has to be modified by particle interactions, 
such as collisions and phenomena caused by space 
charges. In a cylinder in which the magnetic and 
electric fields are parallel to the axis of the tube 
the motion of the particles parallel to the electric 
field remains unaltered by the addition of a constant 
magnetic field applied parallel to the electric field. 
However, the flow of particles perpendicular to the 
Magnetic and electric field lines is inhibited. By 
applying a magnetic field parallel to the tube axis 
the ambipolar diffusion coefficient is reduced by 
(1 + Ton age where Boe is the applied magnetic 
field along the axis and u_ and U, are the electron 
and ion mobilities respectively. 

However, it has been pointed out by Bohm, Burhop, 


(1) 


Massey and Williams that random fluctuations of 
charge density and plasma oscillations may produce 
electric fields which in their turn give rise to drift 


motions across magnetic field lines. This drift motion 


provides an additional mechanism for ionized matter to 
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move across a magnetic field and may, when it dominates 
over collision diffusion, introduce considerable diffi- 
culties into the physics of a plasma in a magnetic 
field. Bohm and collaborators performed experiments 
with a plasma in a magnetic field and came to the con- 
clusion that the diffusion was not consistent with 
collision phenomena. The diffusion across the lines 
was much greater than predicted by collision theory 

and hence an additional mechanism must be present. 


mennents 


studied the diffusion processes in 
a helium discharge in a longitudinal magnetic field, 
by using a long discharge tube. 

The results of Lehnert's experiments showed that 
the radial diffusion agreed with collision theory up 
to a certain critical value of the magnetic field. 
Above this critical magnetic field an instability was 
found to occur and the results indicated a diffusion 
rate much greater than that given by collision theory. 
This gave support of the additional mechanism which 
Supported the enhanced diffusion as stated by Bohm 
et al. 

(3) 


Kadomtsey and Nedospasov showed how this 
effect discovered by Lehnert could be explained by the 
fact that a screw type disturbance will grow in the 


presence of a large enough magnetic field due to the 
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interaction between the current and the magnetic field. 
In the absence of a magnetic field any disturbance of 
the plasma will be rapidly terminated by the increased 
diffusion of particles from regions of excess density. 
In particular, the twisting and wriggling of the current 
channel causes an increase in the flux of particles 
striking the wall at those places where the channel 
comes close to it and a decrease in the flux on the 
reverse side. This has the effect of restoring the 
Original state of the plasma. In the presence of a 
Magnetic field there is an additional force, the inter- 
action of current and magnetic field, acting on the 
plasma. With a helical distortion of the plasma this 
force acts either towards the axis of the discharge or 
towards the walls according to whether the helix is 
right or left handed. When the force acts in the 
direction of the walls there is a build-up of the 
initial disturbance. As the magnetic field is increased 
this force becomes greater and so does the transverse 
diffusion, and along with it, the stabilized action of 
the walls is reduced. Thus for sufficiently large 
Magnetic fields the plasma column becomes unstable. 
Using this explanation Kadomtsev and Nedospasov were 
able to theoretically derive a relation between the 


longitudinal electric and magnetic fields which gave 
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a diffusion coefficient which agreed with the results 
obtained by Lehnert up to and beyond the critical 


Magnetic field. 


THE INSTABILITY MECHANISM 


A physical model of the growth mechanism of the 
helical wave was developed by Hoh and Lehnert ‘*) , 
The argument is begun by assuming a small, quasi- 
neutral, screw-shaped perturbation of positive and 
negative carriers to be superimposed on the steady 


state unperturbed distribution. The general form for 


the perturbation in cylindrical coordinates would be 
F(r) exp[i(md + kz - wt)] 


with F(r) determined by the boundary conditions at the 
wall, >? being the azimuthal angles, k the wave number, 
w the frequency of the oscillation. The perturbation 
is shown schematically in Fig. (la) for m and k having 
opposite signs (a right-handed helix). The applied 
longitudinal electric field E oz tends to lift up the 
ion screw relative to the electron screw. This separa- 
tion is equivalent to a rotation of one screw relative 
to the other, as shown in Fig. (la) and in the cross- 
sectional view [Fig. (lb)]. The charge separation 


resulting from the screw rotation creates an azimuthal 
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Fig. (1) 

Right handed helical perturbation 
for |m=1|. (a) Perturbed-electron 
distribution indicated by dashed 
lines and perturbed-ion distribu- 
tion by solid lines. (b) Cross- 
sectional view showing charge 
separation and resulting perturbed 
electric field E,g. Directions of 
carrier flow due to the combined 
effect of E and B are indicated 


by vectors lian mg® 
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1 
electric field By° The azimuthal electric field 


a 


together with the longitudinal magnetic field Boe 
act on the unperturbed distribution to produce a 
radial flow of particles. The direction of the 
radial flow of particles is determined by the vector 
product r.(E xB) = EyBsz* As can be seen from 
figure (lb) the perturbed azimuthal electric field 
causes the particles to drift to the wall at the peak 
perturbed density and pulls them to the center at the 
minimum perturbed density, thus enhancing the pertur- 
bation. 

On the other hand, the ambipolar electric field 
Eor which is positive will act together with the longi- 
tudinal magnetic field to produce a radial flow of 
particles in the negative azimuthal direction. The 
resulting ambipolar electric field will therefore 
have the effect of reducing the charge separation 


created by the axial electric field E, , Making the 


Zz 
column more stable. 

The effect of a self induced azimuthal magnetic 
field eye) on the column will be to reduce the out- 
ward radial flow of particles caused by the azimuthal 
electric field B, - This is because the induced 


azimuthal magnetic field will act with the applied 


axial electric field to cause a flow of particles in 
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towards the axis of the cylinder. The final flow of 
particles in the radial direction is given by the 


term 
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Instability occurs in the column when the charge 
separation created by the axial electric field is 
great enough to overcome the dissipative effects caused 
by the ambipolar electric field and the recombination 
of particles inside the bulk. It may therefore be 
expected that for a given value of the axial electric 
field E oz there will be a threshold value of the axial 
Magnetic field, But critical, above which the growth of 


the helix is observed. 
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CHAPTER II 


HISTORICAL SURVEY 


IIl.1 -Description.of the isystem 


The system to be studied consists of an electron- 
ion plasma situated in parallel electric and magnetic 
fields. The plasma is generated in a cylindrical dis- 
charge tube and we assume that the plasma is ionized 
weakly enough such that we have to only take into 
account collisions between charged particles and neu- 
trals. We assume that our plasma obeys the perfect- 
gas law, and that although the density of the plasma 
May vary throughout the tube, the density of positive 
and negative particles is the same at every point. 

We assume that internal electric fields may exist in 
the plasma although charge neutrality is assumed. This 
assumption of quasineutrality is reasonable since only 
small differences in the densities can produce very 
strong internal electric fields. 

The apparatus consists of a long cylindrical 
discharge tube and the length of the tube is much 
greater than the diameter of the inner radius, the 
inner radius of the tube being of the order of 1 cm 
while the tube length is greater than 4 meters, so 
that end effects can be neglected. The electric field 


is produced by an electric discharge in the tube and 
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the magnetic field is produced by a magnetic coil. 
The longitudinal electric field is measured with an 
electrostatic voltmeter connected between two probes 
about one meter apart in the positive column of the 


tube. 


2 pie, A: Instability 


The helical plasma density wave (instabilities) 
was first proposed by Kadomtsev and Nedospasov as an 
explanation for the anomalous diffusion observed 
experimentally by Lehnert for a plasma in a longi- 
tudinal magnetic field. In solving the problem 
Kadomtsev and Nedospasov started from two basic 
setsof equations. The first set of equations is the 
equation of continuity for ions and electrons. The 
second set of equations is an equation of motion for 
ions and electrons. They assume that the frequency 
of the collisions between the ions and the neutral gas 
molecules is much greater than the collisions between 
the electrons and the neutral gas and also much greater 
than the cyclotron frequency and the frequency of the 
oscillations of the disturbance. Therefore in the 
equation of motion for the positive ions they are able 
to neglect the effect of the magnetic field. 


The equation of continuity is 
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where "-" refers to electrons and "+" refers to ions. 
iPS nv_ is the particle flux vector, a being tne 
+ + + 


average velocity of either component of the plasma, 
and € is the number of ion-electron pairs created per 
second per electron. The equation of motion which 


relates the flux to the forces is for electrons 


> > 54 > 
fee Deyn + pent. x28 =0... (tae) 


For the ions the corresponding equation is 


> 
i= nk = .00 (ioe) 


til 


insequacilons (ii. 2) tind) (11.3)), D.) is the electron 


diffusion coefficient, ue the electron and ion mobili- 
ties and B and E are the magnetic and electric fields 
respectively. 

Kadomtsev and Nedospasov then solve equations 
CDI G2 prmeandct lh. 3petor e and Te explicitly, take their 
divergences and combine the results with equation 
(II.1) to obtain the two following equations in which 
= 


> 
T_ and T, have been eliminated. 


For the electrons: 
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For the ions: 


on ; Ze a 
rer u,V. (nE) + Ne he Sen0 (eT: 5:) 
where 
' a ' D> 
_ = or, sSO€DW d= —— > - (1156) 
+ Lo + noe eg ucBe 
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In obtaining equations (II.4) and (I1I.5) Kadomtsev 
and Nedospasov used the conditions that V.B = 0. They 
also assumed that the internal magnetic field is negli- 
gible relative to the large uniform longitudinal magnetic 
field which is being externally applied to the column 


and therefore vxB = -0.. 


a) Steady-state solution 


By using equations (II.4) and (II.5) they solve 
for the case of an unperturbed cylindrical plasma in 
a uniform longitudinal magnetic field, where there are 
no gradients of density or electric fields in the axial 
and azimuthal directions. 


The solution they obtain for the density is 
n, (r) = NOG (BSF) P Cites] 


where Bo = a/Rr where OG is the first root of the zero 
order Bessel function Jor R is the radius of the tube, 


Ns is the density along the axis of the tube. 


ve. 


hi 


ake 


|} 
4 


write 


i. an a a 
ay it a) bai 


dolassihall (21a). a ‘Werdbivans ori “ks 
Pare es ea oe re 


svloz yet (2,02) bas (8.2) anoitseps eniew va >: 
at smeskq {so bsbnilys seh lana AS “Re 'pes0- “a 
. 63s sxstitvsweiw ,bfsi2 asi . | 

isixe sdj ai ebiot? 


i* 
a i ; 
al’ 


+ 

Loe : om 7 y war: | = . g hake Dy ¥ oN ¢ a his _ 

Pika ‘A - a 7 fe tat - i 
a, eee 7, 4 ( ee ae a ty eo. ae | 

wy. Pn py) 2) - re, 41+ , he’ ‘Bidet i. a ~ e i 

| yodus sae ae anf : 5X - and bof =) 9 » i, achat ain : a R i, —s +. - a 
ee eas Pel Oe eae 1a iat hixtaalbg 
- a ~ 1 ae ee 

1) ) Mee i, pie Bey cote ata) lise la at or : 
: tess are) i att oe fen . 0, ae ‘ : : 


¥ eal a. ee ‘ a 
a7) ; - yy iy ' ehh x a roe i 
i ‘ 7 , ‘4. we ' mu ; 
oa r | ; bbs 
ae ehh. a Pa pee 
: \ iA “4 > IE ; oh) Y 
7 1 a ; ee eee 
‘ é. : 
7 7 ° » 


LS 


They then determine a value for the ambipolar 
electric field by requiring that the radial flux of 


ions is equal to the radial flux of electrons 


Ic tae . (L163) 


The value they obtain for the radial electric field is 
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Now 
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Te = Ge BJ (8GY) F (LE« LO) 


The radial electric field becomes 


Da J, (Bor) 
E = 
OF eee ee fe) J (Bor) 
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The radial electric field is positive, that is it points 
from the cylinder axis towards the walls of the tube. 
This is what you would expect of the ambipolar electric 
field. Also since J, (Bor) is zero at r = 0, along 
cylinder axis, the value of the radial electric field 

is also zero on the tube axis. For symmetry reasons 


this is what you expect. 
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b) Time-dependent solution 


The authors then solve for the conditions under 
which the plasma becomes unstable, that is the condi- 
tions under which the anomalous diffusion as observed 
by Lehnert will be dominant in the tube. To solve the 
problem the authors have chosen to consider a small 
perturbation in their calculated steady state values 
for the density and space charge potential. The form 


of the perturbation they choose is 


Ber) exp l[i(me a kz) = wt) 7) 2, CT raleZ ) 


where @ is the azimuthal angle, k is the wave number, 
w is the frequency of the oscillation. Both @¢ and k 
are real numbers while w is a complex number. 

The density and electric fields can be written 


as 


Gileei >) 


where No and E, are the steady state values of density 
and electric field, ny and E, are the perturbed values 
of the density and electric fields. Equations (II.13) 
are then substituted for n and E in equations (II.4) 


and (II.5) and only terms of the lowest order in the 
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perturbation are kept. Equations (II.4) and (II.5) 
become: 
For electrons 
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Equations (II.14) and (II.15) contain only the perturba- 
tion terms and do not contain any steady state terms 
which have been grouped together and equated equal to 
zero separately. 


nj and a are written explicitly as 


n, = f(r) exp[i(kz + md - wt)] 
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Since we demand that the density be zero at the tube 
walls, f(r=R) = 0. 


They define a new variable 


g(r) = : (II.18) 


Since g(R), which is essentially the amplitude of the 
potential of the electric field at the boundary, must 
be finite and since n(R) has finite values, due to the 
boundary conditions at the walls, 2%(R) must then be 
finite. Using condition (II.18) they substitute 
equations (II.16) into equations (II.14) and (II.15) 


obtaining: 


For electrons 
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For ions 
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The authors then assume that the radial dependence of 


the perturbation has the following form: 


| 
oe J, (BOY) 
COPS 22.) 
| 
where By = (a,/R), where Oy is the first root of 


the first order Bessel function, n' and V' are cons- 
tants independent of the radial positions. Each of 
equations (11.19) and (I1I.20) contains a real and 
imaginary part which can be separately equated equal 

to zero. Equations (I1I.19) and (II.20) contain five pa- 
rameters, the wave number k, the frequency of the 


oscillation w, the azimuthal mode m, the applied 
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electric field Eos and the applied magnetic field Bas 
Using the conditions that the imaginary part of the 
frequency must be greater than zero and that the 
derivative of the frequency with respect to the wave 
number must also be zero, they are able to derive a 
relation between the critical electric and critical 
Magnetic field for the m = 1 mode. 

The reason for demanding that the imaginary part 
of the frequency be equal to zero can be shown by the 
following calculation. The condition that the column 
becomes unstable is that the helix grows with time. 
Therefore if we look at the time dependent portion of 
the instability the real part must be greater than l 
for the helix to grow. The time dependent portion of 


the helical instability is given by 


exp[-iwt] . | CLE 7233) 


Now 


Oo = Of + LHW. : Cite 4) 


Putting (I1I.24) into (I1I.23) we get 


exp [-iw t] exp[w,t] 5 CErS25)) 


We can see from relation (I1I.25) that the condition 
that the wave grows is that expl[w,t] is positive. 


Therefore Ws must be greater than zero for the helix 


oi Yd aworla oil aso O08 of Lsupe od ROUpST3 | 
amuioo add ted3 avittbado oat : :Setvoten pat ” E 


i asda able ad seuMn ilies ‘a sds yoel | 
to wichita agsbasqeb smid sit .woxp ‘oo xilod afi” | 
yd ata ei yiilidstans Leo tied da i‘ 


(ES. 1D) | aay bites 
| crea -_ 


top ae IZ) ojni- pte 
ny 


(avery Le 7" 


i 4 Charente } 
‘ a i 7 ie 


to grow. If w;= 0 this will define the boundary 
between the region for which the helix grows with 
time and decays with time and hence the boundary 
between the stable and unstable regions. All points 
on this boundary would be critical points. 

Their calculated theoretical critical points 
agreed with the corresponding experimental critical 


points obtained by Lehnert ‘2? , 


(5) 


Johnson and Jerde extended the theory of 
Kadomtsev and Nedospasov in that the diffusion and 
Magnetic field interaction terms are included in the 
ion equation of motion, and also in that the equations 
are solved in a more rigorous manner without any prior 
assumptions regarding the form of the radial depen- 
dence of the perturbations of density and potential. 


The equation of motion for the ions used by 


Johnson and Gerde is 


> > + a 
gt efi mBlde Davna- (1, B) p, - (TE26) 


Equation (I1I.26) differs from equation (II.3) in that 


it includes the diffusion term for the ions D,Vn and 


> > 
the magnetic field interaction term Hens <B)i. 


Following the same procedure as Kadomtsev and 
Nedospasov, Johnson and Jerde used equation (II.26) 
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the equation <ofrcontinuity for dons: equation “(II.1). 


The equation obtained is 


on J = re? : Eg eo 
arm tae c= u,V. (nE) + DV Dot Hu, V. (n (B x E) ) 
Pe ie So are 
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Equation (11.27) ditters from equation 3.11.5). im that 
it includes the diffusion term and ahs terms due to 
the magnetic field interaction. 

Eliminating the radial electric field between 
equations (11.27) and (11.4) they obtain an equation 
for the density of the unperturbed plasma. The 


equation they obtain is 


= (x =) + Ben = uO! (IT.28) 


The solution to this well known equation is a Bessel 


function of order zero 


n= NoJG (BGY) P ele) ) 
where 
Cech aie 
g¢ Be baie ies : (irs30) 
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and No is the density along the cylinder axis. 
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The only difference between equation (I1I.29) and 
equation (1157) is the definition of the first zero 
Bo Thaseoccurs, asia result: of the, inclusion of the 
diffusion and magnetic field terms into the theory 
for ions and results in new coefficient for ny being 


formed when Ey is eliminated between equations (I1.27) 


z 
and wiChi. 4) < 
They then .calculated the radial electric field 


using condition (I1.8). The value obtained is 


be (EL a31) 


Using equation (I1I.29) this can be rewritten as 


D -D J, (2.45) 
E = st ie Soni 
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or Tie ei fe) J, (8Y) 


The radial electric field (II.32) differs from that of 
Kadomtsev and Nedospasov only through the mobility 

and diffusion coefficient terms and in the value of 

the first zero of the zero order Bessel function. It 
can be seen by comparing (11.25) and (II.11) that the 
diffusion coefficient for the ions is now present whereas 
it was not present in equation (II.11). 


Equation (II.25) can be rewritten as 
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It can be seen by looking at equation (II.33) that Kor 


is positive provided the relation 
“Les CID. 34) 


holds true. 
However, it would seem that for large values of 
Bos the relation (11.34) would not be true, therefore 


the radial electric field would be negative and point 


in towards the cylinder axis. This observation was not 


pointed out by Johnson and Jerde. 

Johnson and Jerde solved for the condition for 
the onset of the helical instability using the same 
perturbation method as Kadomtsev and Nedospasov. The 
inclusion of the diffusion and magnetic field terms in 
the equation of motion for the ions resulted in a dif- 
ferent equation of continuity for the ions. That is 


in replace of equation (iI.20) they used 
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However, Johnson and Jerde did not make any 
assumptions regarding the radial form of the per- 
turbed density n, (r) and perturbed potential V, (4) 
as did Kadomtsev and Nedospasov. Instead they made 
Hankel transformation of the complete equations of 
continuity of the perturbed ions and electrons, 
equations (I1I.35) and (I1I.19) respectively. The 
independent variable r was transformed into Oss 
where Jyq (52/R) = 0, where R is the tube radius. 

n, (r) and V, (4) were transformed into the functions 
Pe) and Litas) They then were able to eliminate 
Fo.) from the system of equations and the remaining 
system of algebraic equations could be written as 


L(a.) + } 


Eaton) = 0.4; (EL 2354) 
i=) 74 J 


where the matrix elements aye are known functions of 
axial magnetic field and constants. Due to the extreme 
complexity of ee only the first term in the matrix 


was used 


Ne Se = 0. 0) 


Equations (II.35b) were solved, using the same conditions 
as Kadomtsev and Nedospasov, for the critical electric 
field in terms of the critical magnetic field. They 


then put these points ona graph (fig. 2). One can see 
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Fig. (2). Stability criterion for the helical 
mode in helium as a function of the 
longitudinal magnetic field. 
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by looking at the resulting curve that the larger the 
applied axial electric field the smaller is the re- 
quired axial magnetic field necessary to make the 
instability grow. This is consistent with the theory 
for the mechanism of the instability as pointed out 
earlier. 


(6) 


Paranjape et al. took the same approach as 
Johnson and Jerde, but modified the problem by consi- 
dering the effect of a small self induced azimuthal 


Magnetic field. The azimuthal magnetic field is simply 


related to the electric current by 
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where ve are the z components of the drift velocities 
which they assume to be independent of the radial and 
axial directions. They also defined an average plasma 


density given by 


r 
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Therefore the azimuthal magnetic field can be approxi- 


mately written as 
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Following the same procedure as Johnson and Jjerde they 
obtained an equation for the steady state plasma. The 


equation they obtained is: 
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B can be rewritten as 
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Therefore B can be written as 


where g is defined by Johnson and gerde. ‘>? ” aie & | 
equation (II.40) we put the azimuthal magnetic field 


term equal to zero (A = 0), the solution we obtain is 
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No is the density at the origin. This can be re- 


written as 
n = noJ, (BY) , 


which is the same solution as Johnson and Jerde: this 
is what you would expect to obtain. 

Using the method of perturbation the next appro- 
xXimation to the solution of equation (II.40) they 


obtain is 


N = NG LI, (YP) oe a,5,(y,0)1 F CEE) 


where pas 5.52 is the second zero of the Bessel func- 
tion of order zero and ay is a small number. The term 


a Jo fy /) increases the density at the center from 


1 
NOG (Yo?) and decreases the density near the boundary 
wall. 

The radial component of the internal electric 
field was determined by requiring that the radial 


flux of electrons is equal to the radial flux of ions. 


The resulting radial electric field they obtained is 
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It can be seen from equation (III.48) that if we 
put the azimuthal magnetic field equal to zero B, = 0 


o 


they obtain the same value for the radial electric 
field as Johnson and Jerde. It can also been seen from 
equation (III.48) the effect of Ss is to reduce the 
radial electric field. This is what one would expect 
Since the induced magnetic field increases the density 
towards the center of the bulk and reduces the density 
for the larger radial values. Therefore the flow of 
electrons to the wall is decreased and hence a smaller 
radial electric field is required to keep the flow of 
ions to the wall equal to the flow of electrons. 

Using the same method as Kadomtsev and Nedospasov 
they derived a relation between the critical electric 
and critical magnetic fields for various strengths of 
the self induced azimuthal magnetic field. They then 
plotted curves of the critical electric field against 
the critical magnetic field for various strengths of 
the =induced magneticrrreld (fig. (3)). It can be seen 
from the graph that for increasing values of the 
induced magnetic field increasing values of the applied 


axial magnetic field is required to cause the instability 
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to grow for a given critical electric field. Therefore 
the effect of the self induced azimuthal magnetic field 


is to make the plasma more stable. 


€))j me loinstabilitys Fanite Amplitude 
(7) 


Holter and Johnson have expanded the theory 
as developed by Johnson and Jerde by superimposing a 
finite amplitude helix upon the steady state density 
and potential instead of a perturbation of the steady 
state density. They first calculate the functional 
form for the density of the background plasma including 
the effect of the finite amplitude helix on the back- 
ground plasma. They then determine the stability con- 
ditions for the plasma. 

When solving for the background density, steady 


state case, they assume that the density and space 


charge potential have the following form 
n= Re {Noho (r) + N,f(r)exp[i(wt + kz + mo) ]}. 


Ay is the background plasma density which is affected 
by the amplitude of the helix. 
V = Re Lp. (x) + V,g(r)exp[i(wt + kz + mo Foye bh, (LLL 49) 


where Ne is the density along the cylindrical axis, Ny 


is the finite amplitude of the helix. ee is the space 
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charge potential of the background plasma which is 
affected by the finite amplitude of the helix, Vi is 
the amplitude of the helical potential. 6 is the 
phase shift between the helical density and the 
helical potential and is a real number. Equations 
(II.49) are then substituted into the equation of 
continuity for the electrons equation (II.4) and for 
ions equation (11.27), keeping all terms including the 
non-linear terms which are terms with coefficient pro- 


ducts NjVj- The resulting equation obtained is 
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Equation (11.50) differs from the corresponding 
equation as derived by Johnson and Jerde in that they 
include the generation of higher harmonics from the 
first harmonic term. This can be seen by looking at 


the terms which multiply 


Expiltege (it) tekZeatam®) fF 6) N,V, . 


It can also be seen from looking at equation (II.50) 
that the inclusion of these higher order terms affects 
the zero harmonic equation. This can be seen by look- 
ing at the two terms which multiply NjV)- It should 
be pointed out that the non-linear terms or second 
harmonic terms in equations (II.50) are incomplete 
Since they do not include the higher harmonic terms 

in the definition of their density and space charge 


potential equations (11.49). This omission implies 
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that they do not incorporate in their theory the genera- 


tions of higher harmonics by the finite amplitude first 


harmonic. 

To obtain the steady state solution they look 
at the zero harmonic portion of the equation setting 
it equal to zero and assuming that wit << 0. The zero 


harmonic portion of the equation becomes assuming m=1 
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To solve for the density AG they first eliminate Ey 
between the ion and electron equation in (III.51) and 


obtain the following equation 


2 ee fg) = 
ee nace ee) — 0 (Tre 52) 
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Equation (11.52) differs from the corresponding equa- 
tion in the theory of Johnson and Jerde, equation 


(II.28), in that it contains the term 


2 fg 
OV png 


r 
One can see immediately by looking at equation (II.53) 
that this term arises from the inclusion of the non- 

linear terms in their equation (I1I.50). This term was 
neglected by Johnson and Jerde since they considered 

their helical density distribution to be perturbation. 
If you put ¢ = 0 you obtain equation (11.28) which is 


what one would expect. 
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Holter and Johnson solved equation (II.52) by 
writing it in integral form and using certain rela- 
tions between Bessel functions. They also used the 
boundary conditions that the density be one at the 
origin (h. (0) = 1) and zero at the tube wall GS: Nie 0. 


The equation they obtained for the density Ae is 
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Equation (I1.54) is a non-linear integral equa- 
tion for ha and was solved by Holter and Johnson 
using a digital computer. It can be seen by looking 
at equation (11.54), if » is put equal to zero, which 
is the same a neglecting non-linear terms, we get the 
well known density solution which is a zero order Bessel 
function: 

Holter and Johnson plotted a graph of relative 
density against radial location for various positive 
values of » (fig. (4)). It can be seen from their 
graph that for ¢ equal to zero the density profile is 
that of a zero order Bessel function. As $ becomes 


very large, the density profile is given approximately 
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Frg. 


Radial Location 5 
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(4). Radial density profile of the 


positive column for 9=0, 3.27 


and 8.92. 
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by a first order Bessel function (i.e. J, (8,¥)) with 
a deviation occurring near the origin. 

Holter and Johnson then calculate the radial 
electric field not by requiring that the radial flux 
of electrons and ions are egual as Johnson and 
Jerde did. Instead they eliminated the — term from 
equation (II.51). They assume the following form for 
the radial portion of helix 


J, (8,1) 


£ (xr) =J, (61) g (xr) = h (rt) Z (2 235) 
Oo 


The resulting equation they obtain for the radial 


electric field is 
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dh a6 Jj (8,4) 
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V are the terminal energies for the ions and electrons. 


ire 
k,T_ 
Vv = —+ 
i e , where Ms are the ion and electron tem- 
peratures. 
Looking at equation (II.56) it can be seen that 
if we put ¢ = 0 the solution for the radial electric 


field is the same as that obtained by Johnson and Jerde. 
Holter and Johnson using the first harmonic por- 
tion of equation (11.50) and following much the same 
procedure as Johnson and Jerde obtained a relation 
between the critical electric and critical magnetic 
fields in terms of the parameter $¢. It can be seen 
From Eig. .(5) "esata for a‘given critical electric field 
the larger the value of $» (that is the larger the ampli- 
of| the helix) the largemis=the critical ‘magnetic field 


required to make the plasma unstable. 


ETE 3a em 2 instability 
(8) 


Holter and Johnson extended their own theory 
by investigating the conditions whereby the helical 
instability (m= 1 mode) could excite higher modes in 
the plasma column. They investigated under what con- 


ditions on m = 2 mode can be excited in the presence 


of a finite amplitude helical mode (m = 1). They 


pate 09k ed 60! ae ean POE tian a : a: 
pixtnste isibex pe rol nogaulos: afd =o dug 9 Mass: 
-sbiet bas noendod ve oankesdo ‘g46dt 86 ‘amee eae ek bias 


¢ > 7 
a 


omsd: ect doom priwollod bus (tee ) eimies. 26.) 
sobbed 6 bénissdo sbast. bas aoendiot a6 
DitonpBA Iepitivo Sas oizdesie fssadizo orit 
mese od S169 arb aianew eS ons to aries ries able 
bleit ckxtoels Isoidixo gevip 6 X02 art @) «pid nt 
-ilqms odd xopisit st ei ted3) 4 to ouley ort oaneke ce 
Siok? oissapam Lsoisixzo sdt ek “zepiel! ‘eds (xifteod. ee 
-oidstenu smesigq- dt asm ot . 


a ala Awo xiorls sonata 
lepkted ‘eds ydousdo 2 BE 
ak aebomn xedpid — — 

= fOD eit xebi bases 


ele gre Uc) 


—- sUns table 


107° 107.03 10¢ 3 10° 3 10210 10° 107 10 


2 
El ha hagegpa 


Stability criterion for the finite m = l 
instability as a function of the longitu- 
dinal magnetic field for different values 
of the amplitude factor ®6. 
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consider the m = 2 mode to be a perturbation of the 
m= 1 mode. They express the density n and space 


charge potential V in the form 
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As in the previous theories the axial wave number i 
and the azimuthal wave number m and the phase angle 
On are real, while the frequency On is complex. “The 
amplitudes Noe VES Ny and Vi are finite and N, and V5 
are considered small quantities such that products 
NoV, are negligible. 

They calculated the background density profile 
and radial electric fields and obtained the same values 
as they obtained previously. This is reasonable since 
they consider the terms NoV> are considered small and 
hence the m = 2 mode will not affect the background 
density or space charge potential in the lowest order. 
They then derive stability conditions for the m = 2 


mode, and compare their calculated critical fields with 


experimental observations. 
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Daugherty and Ventrice 


(8) 


extending the theory 
of Holter and Johnson determined the radial form 
Of thenparticle density, bys consadering ithe effect of 
a finite m = 2 mode on the previously calculated 
background density, which had used a finite m= 1 
mode. They use the same form for the density n and 
the space charge potential V as in equation (II.57) 
except that they consider the term No and V5 to be 


finite. They also assume the radial form of the helix 


to be 


Peat) = J (Bt) 
Ciro 6) 
Jig (Bnd? 
cas) = as ’ 


where hy (r) is the steady state density and will be 
affected by the inclusion of a finite m = 2 mode. 
Substituting equations (II.57) into the equation of 
continuity for electrons, equation*® (11.4), and for 
ions equation (11.27), they obtain an equation which 


corresponds to equation (11.50). This equation 


differs from equation (II.50) in that it contains terms 


of order N,V. The second harmonic equation is com- 


plete and also partial third and fourth harmonic terms 


are included. Following the same procedure as Holter 
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and Johnson they eliminate the radial electric 
field from the zero harmonic equation and the result- 


ing equation obtained is 


O 
N,V 
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The only difference between equation (II.59) and the 
corresponding equation in the theory of Holter and 


(7) 


Johnson , equation (11.52), is the additional term 


ea 12682) ; (II.61) 
O 

Looking, até the sdefinition of 5 froma (lEo6O0rnpt tins 
clear that expression (II.61) arises as a result of 
the inclusion of the finite m = 2 mode in the density 
and space charge potential. This is not included in 
the theory of Holter and Johnson. So by putting d5= 0 
we get the identical equation as (11.52). It should 
be pointed out that ? = ¢ when m is put equal to l 
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By requiring that h, (0) = 1] and he (R) ae etic 
authors were able to express equation (II.59) in the 


integral form as 
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where x_ = Bur: 


Looking at equation (II.62) one can see that if 
you put b> = 0 which is neglecting the N5V5 terms, or 
considering the m = 2 mode as small, one gets the 


(7) 


solution as obtained by Holter and Johnson equation 
(1t.o4) 8 LL you puG-bDotn oy and b> equal to zero, that 
is to consider the m = 1 mode a perturbation of the 
background density, you get the solution of Johnson 

and Jerde equation (10.29). Finally, if you put dy 
equal to zero and consider o> as being finite you 


obtain the case of a finite m = 2 mode occurring with- 


out the presence of a m=l1 mode. 


42 


ae pie | 


s - 


oo * hh ane 


($8.22) 


ti tet ssa nso ante to atilp sepia nuigintiaeattcala! OE | 
TO .snaed Veh oft paismelpen wi siokdw 0 = 66 dug voy a 
| ae ae, eo Linh ae bao 5S a eee : 
rottevpe ") noanriot bas xedloH yd benksddo a5 


sais ores oF Laupe of "Gna pp tisod tuq Moy XE _(beat} a 
| nat bela «eer nee 


43 


Equation (II.62) was solved with a digital 
computer by Daugherty and Ventrice and the results 
obtained were plotted on four separate graphs show- 
ing the density against the radial position in the 
tube for various values of ¢, and 9,- 

The first graph (fig. (6)) shows the density 
radial profile when 5 is put equal to zero for 
various values of ¢,. The curves obtained are the 
Same as those obtained by Holter and Johrison ‘7? neal 
fig. (4). This is what one would expect since putting 
b> = 0 reduces equation (II.62) to that of equation 
(11 .54)0) sehooking at ithe curvesin fig: (6) we canusee 
that for $, equal to zero the radial form for the 
density is a zero order Bessel function and as oy is 
increased the radial density function approaches a 
first order Bessel function, with a deviation occurring 
near the origin. As oy becomes large there is an off 
axis peak occurring at approximately r/R = .46 with 
no Gip/oceurring near the origin. 

inthe second graph (fig. (7)) “they show tne case 
of oy equal to zero and 5 Varying, Vlhis 1s thea case 
of a pure m = 2 mode occurring without the presence of 
am = 1 mode in the column. They found that as > 
became large the radial density profile approached a 
second order Bessel function with a deviation occurring 


near the origin. For large values of ¢, an OL axis 
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Fig. (6). The particle density vs radial 
position, in the presence of a 
finite m = 1 mode instability. 
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The particle density vs radial position, 
in the presence of a finite m = 2 mode 
instability. 
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peak in the density occurs at approximately r/R = 0.60. 
POL sita li wolues of > it was found that the radial 
function decreased more rapidly than a zero order 
Bessel function. Unlike’ the previous’ graph the profile 
shows an off axis dip before it rises to its maximum 
value. 

The third graph (fig. (8)) shows the case of a 
constant on = 335 /and 5 varying }erom' +6.0 tora 
Maximum value of 22.7. (Looking at 5 in equation 
(II.60) we see that it is possible to have a negative 
5 if the phase angle 65 is chosen such that sin O5 
is negative). It can be seen from the graph that as 
5 is increased positively, the off axis peak increases 
in magnitude and shifts toward the tube wall. As 5 
is increased negatively, the off axis peak decreases 
in magnitude from the 5 equal to zero case and shifts 
toward the column axis. But no dip seems to occur near 
the column axis. 

The last graph (fig. (9)) shows the case of a 
constant 5 fixed at a value of 20 while dy is varied 
from ooo. Goon. "aS oy is increased positively the 
off axis peak increases in magnitude from the o4= 0 
case and shifts towards the column axis. In addition 
a dip near the column axis becomes less pronounced 


and shifts towards the column axis, but is eventually 
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The particle density vs radial position, 
in the presence of both the finite m= l 
and finite m = 2 modes of instabilities, 
for fixed m = 1 mode amplitude factor 

and variable m = 2 mode amplitude factor. 
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Fig.(9). 
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The particle density vs radial position, 

in the presence of both the finite m = l 
mode and finite m = 2 mode instabilities, 
for fixed m = 2 amplitude factor and varia- 
ble m = 1 amplitude factor. 
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removed for large d,- AS $] is increased negatively, 
the off axis peak decreases in magnitude and shifts 
towards the wall. The dip in the density becomes 
more pronounced and also shifts towards the wall of 


the tube. 


5 er lee) Experimental Evidence of a Pure m= 2 Mode 


The case which is of most interest is that of a 
pure m= 2 mode (fig. (7)). It does seem that a 
m = 2 mode could not occur without the presence of a 
fully developed m = 1 mode first. Ventrice and 


(10) 


Massey performed experiments using a very short 
wide tube with a length to radius ratio of 6. This 
differed from previous experiments where a long narrow 
tube had been used with a length to radius ratio 
greater than 400. The results they obtained were 
Similar to those in fig. (7). Under the experimental 
conditions they used it appears that the conditions 
necessary for the development of the m = 1 mode could 
not be met whereas the conditions for the m = 2 mode 
could be satisfied. 

One of the experimental graphs obtained by Ventrice 
and Massey is given in fig. (10). The resulting curves 
do resemble those in fig. (7) giving evidence for a 
fully developed m = 2 mode in the tube without a 


developed m = 1 mode. It can be seen from the graph 
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Fig. (10). 


Radial POsition, = 


The relative particle density (n/n_) 

vs the radial position (r/R) for various 
values of the applied longitudinal mag- 
netic field B with a constant tube 
current of 100% ma. 
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that for larger values of Bae the larger is the 
maximum value of the density curve. This is also 
in agreement with the curve in fig. (7) because as 
5 is increased so is the maximum value of the 
density curve. Looking at equation (II.60) we can 
see that 5 depends linearly on Be sOutchat ).as Boe 


is increased so is o> 


ee) have pointed out in their 


Huchital and Holt 
paper that a transverse component of an externally 
applied magnetic field has a substantial effect 
upon the critical magnetic field strength. They 
even suggest that a large enough transverse com- 
ponent of the magnetic field could make the m = 2 
mode rise in the tube before a m = 1 mode is 
developed. 
It was with this in mind that they then per- 
formed an experiment using an applied axial electric 
field and a magnetic field with a large axial component 
on a small transverse component. It was observed that 
for a transverse magnetic field of 70.5 gauss and an axial 
magnetic field of 1880 gauss that the m = 1 mode 
was less stable than the m = 2 mode, that is to say 
the m = 1 mode grew more rapidly than the m = 2 mode. 
However, for a transverse magnetic field of 81.5 gauss 


and a longitudinal magnetic field of 1840 gauss the 
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m = 2 mode grew before the m = 1 mode. Therefore it 
does seem possible that if the applied magnetic field 
has a large enough component in the transverse direc- 
tion it is possible for the m = 2 mode to grow in the 


tube before and independent of the m = 1 mode. 
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CHAPTER fi 


MODIFIED STEADY STATE SOLUTION 


The purpose of the following calculation is to 
find under what conditions, if any, we can have a 
pure m = 2 mode growing in the tube before a m= 1 


(13) have 


mode is developed. Ando and Glicksman 
pointed out that this situation may be possible if 
the azimuthal magnetic field is taken into considera- 
tion. In our calculation,therefore, we will take 
into consideration the effect of the azimuthal mag- 
netic field. Also we will reduce the effect the 
walls have in causing the instability by considering 
a wide tube, instead of a long narrow tube as done 

by previous authors. We will also consider that most 
of the recombinations and generations occur in the 
bulk. The major factor affecting the radial density 
profile in the steady state solution is the induced 


azimuthal magnetic field instead of the effect of the 


walls as in the case of previous calculations. 


Lit.l Formulation of whe Problem 


We start our calculations from the same two 
basic equations as used by Johnson and Jerde. The 
first equation is a steady state force equation which 


relates the flux to the forces and can be written as 
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"=—""refers to the electrons and “+" refers to the 


ions, where n is the number density of ions or 


> 
electrons, T. are the particle flux, and are defined 


+ 
> 
by i = nv_ , where V_ are the drift velocities. 
+ + + 
De A= kel Je = ny are the diffusion constants, 
+ 4 BE + F 


and Pe are the mobilities. i are the temperatures 
and ve = ae are the thermal Brergtes expressed 
in electron volts. E and B are the actual electric 
and magnetic fields. Both E and B include the inter- 
nal and external fields. The form for the fields in 


the steady state problem are as follows: 
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Foz and Bo, are the externally applied fields, while 


E, (rv) and B,(r) are the internal fields. B,(r) is 


v 
the induced magnetic field in the azimuthal direction 
and is ‘simply related to the electric current by the 


equation 
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+ a 
where ve are the Z components of the’ drift velocities. 


Let us now examine the terms in equation (III.1). 


The first term is 


This is the flux of particles which arise from the 
density gradient and its direction is the same for 
both types of particles, although its magnitude is 
MUCK larger tor electrons jthan, for sons. In the 
steady state case there is only a density gradient 
in the radial direction and therefore the flux of 
particles will be in towards the center. 

The next term is F¥ wank. Toaiseis. them lux oF 
particles due to the electric field. The applied 
electric field (III.2) causes the particles to move 
in opposite directions in the axial direction. 
However, the radial electric field acts in such a 
Manner as to retard the flow of electrons to the wall 
while it enhances the flow of ions to the wall causing 
ambipolar diffusion. 

As > 

The last term is +# He (Te aoe nS es tne £1ux 
of particles which arises as a result of the Lorentz 
force of (V x By where Ve is the velocity, due to 


the current and the magnetic field. This term causes 
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the electrons to gyrate about the Bere lines ina 
right-handed sense while it causes the positive ions 
to rotate in the opposite direction. In the steady 
state case if you take the radial component of this 
expression you obtain 


Fon! Peride “HY cued eBid dl] tie) 
= z + F 9 


where = is (che Sux oLr particles an tne-z Gizrcection. 
Since the electric field is applied parallel to the 
axial magnetic field we will assume that they are both 
applied in the positive axial direction. Therefore 


the ions will move in the positive z direction while 


the electrons will move in the negative z direction. 


Now 
> > 
Das ek et) S 
+ + 
- cs = F nv : 
+ + 


Equation (III.4) can be written as 


Eg Pate ee 
TRS etre ex BIS =e enV BAe (III.4a) 
cat ites + +? 


It can be seen from equation (III.4a) that the effect 
of the self induced azimuthal magnetic field is to 


reduce the outward radial flow of particles and 
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therefore make the column more stable. If instead of 
taking the radial component of the expression we take 
the azimuthal component, we get 


Fhe Demi hab ke i: (III. 4b) 
i ps Ee aenOz 


Now since the outward radial flow of particles is 
ambipolar ra is" positive for both’ dons and electrons: 
Therefore it can be seen that the effect of Bae ES eto 
cause the electrons to gyrate about the Eis Lines iin 
a positive manner (right-handed sense) while it causes 
the ions to rotate in a negative direction. 

Our second basic equation is the equation of 
continuity 


9 (n-n) e he 
Teresi + Vea = J (nen ye, (io) 


where 1/& is the mean bulk lifetime for excess carriers, 
and n is the density at which no net recombinations or 
generations take place. 

Equation (III.5) states that the time change 
in the net number of particles available for recom- 
bination or generation in a small volume plus the 
divergence of the flux of particles is equal to the 
net number of recombinations or regenerations in that 


volume. 
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We next put our value for [T. into equation (III.5) 
+ 


and obtain for our continuity equation 
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We have used cylindrical coordinates for reasons 


> 
of symmetry in solving for [I The z axis is taken 


t° 


along the direction of applied electric and magnetic 
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fields. We have neglected [B, (x) /By 1° in comparison 


with unity. Here B,(r) is the self induced azimuthal 


> 


magnetic field. 


EEY.2 ‘steady State Solution 


In the steady state solution there are no varia- 


oy) 


tions in) the: density with time. .Therefore 9(n-n) /9t=0. 


In view of the symmetry and our choice of cylindrical 
coordinate the density distribution will be in the 


radial direction, therefore we need only consider 


ourselves with the radial component of equation (III.8). 


Therefore using our steady state condition and 


remembering that V.B=0, only the r dependent terms from 


equation (III.8) remain 
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Using equations (III.9) we eliminate the radial 
electric field E (©) and obtain one equation for n(r) 


in terms of the induced magnetic field B,(r) and other 
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constant terms. 
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Using definitions (III.7) and a u_kT_/fe and 
+ + 
dividing tthrough by? (me + me equation (III.10) 
becomes 
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Multiplying both sides of equation (III.11) by rdr 
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Defining the following dimensionless quantities 


x= q(r) = 


222 
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where R is the radius of the tube and n is the aver- 
age density of the plasma. Using equations (III.14) 


equation (III.11) can be rewritten as 


x pet x 
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Equations (lI Iai )@asmeatnone lamearl equation and 
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solved by penpene eee when a4= Ofeewthes solution 
obtained is 
2 
Tp (x) 5 ' COI Esse: 
eve 
Lao ee) ] 


e8 asttitwer sd’ ‘ie (UL. 15t) notteupe al 


oe 
; ‘ 
Me 7 


(@£.472) shan = (2)e)) be 1 sie] Bi. ‘. 


where ZnO. =Eo“7T gi and x/X * 1,fand Q, is the 
Deeasmaidensity at the origin. In putting 4 equal 
to zero Bennett has neglected the recombinations 
and generations which occur in the bulk of the 
plasma. We attempt to include the effect of re- 
combinations and generations in an approximate 
manner. 


Let us assume that the term 


x 
al 
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0 
in equation (III.15) is small. This is reasonable if 
Oy is small. We will consider the solution to equa- 
tron (hiss) tov besas perturbation or the Bennetr 
solution, (equation (III.16)). We write the solution 


as 


q(x) = Ap (X) ony ela aer (i fae 7) 


where*y (x) is the ~perturbed* terms "We Let°y(x) be 
proportional to Oy and expand the coefficient of ot 
in a power series of x, keeping only the first five 


terms. Therefore y(x) can be written as 
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where the coefficients Agr A Ayr Ay have to be 


ioe 
determined. 

To determine the coefficients in the expansion 
fOr’ y¥(x)Swewsubstiltuge  (#ii.17) intor (LI1I.15) using 
(III.16) and (II1I.18), keeping only the lowest order 


terms in 4 (i.e. neglecting terms of order Gain 


Equation (III.15) becomes 
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Using equation (III.19) and equating the co- 
efficients of like powers of x equal to zero, we can 
solve for the coefficients in the expansion for (x). 
In determining the coefficients we have used the 
binominal expansion, since we have assumed (x/X.) < iy 
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It follows that 
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The coefficient Ag cannot be determined from equation 


(III.19), but will be determined below. 


Ties) Catculration of AG 


To calculate Ay we write down the equation which 
states that the total number of generations inside the 
plasma is equal to the total number of recombinations 
at the walls. However, Since we are considering a 
wide tube, as pointed out earlier, the effects of the 
walls on the density in the bulk will be small. The 
equation is 

R 

| E(n(r)-n)rdr = RS(n(r)-n) , Cries 

0 
where S is the surface recombination velocity and n 


is the density at the wall at which no net recombina- 


tions or generations take place. 
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It should be pointed out that although equation 


(ETE 22) 1Ssean, exact, equationgour solution for a AG 
will be crude. The reason for this is that we keep 
terms up to powers of x4 in ous uadensa@ty* distribution: 
Therefore a truly exact value for AAG HS enol. ob— 
tained. However, our calculated value will give us 
an approximate correction to the Bennett solution on 
and near the axis. 

It can be shown by physical arguments that a A, 
must be negative. The generations inside the bulk 
will occur when the density of the plasma is less than 
n. Now since the Bennett solution is a maximum at the 
axis and a minimum at the walls, and since we insist 
on a net number of generations in the bulk, these 
generations will have to occur away from the tube 
axis. However, the Bennett plasma density along the 
axis, Qo: is greater than n therefore along the 
axis there will be net recombinations and the plasma 
density will decrease from the Bennett value, Qo: 
Therefore we would expect from physical arguments that 
the first term in our density perturbation to be nega- 
tive. 

Rewriting equation (III.22) in terms of variable 
x and aoing the required integration keeping only terms 


up to a,A,, we obtain a very approximate equation 


which can be solved for OA, 
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Se AGF aA, + aA, iu] (2EEiv23) 


where 


In obtaining equation (III.23) we have used the fact 


1 = A3 = 0. 


We rewrite equation (III.23) expanding the term 


that A 


[l + (1/x,)71 7 in a binominal expansion keeping 
terms up to (xe This expansion is allowed since 
1/x,, <1. Also, we make substitutions into equation 
(III.23) for the coefficients a A5, a Ay using equa- 
CLOUS MCP) . 


Equation (III.23) becomes 
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We will keep only the dominant terms in equation 
(III.24). Since we are considering the effect of the 
walls to be small we can consider Pp to be small, such 
that we can neglect terms of order op, play | relative 
to unity. We also neglect Jo, | (1) relative to unity. 


Therefore equation (III.24) can be written as 


Ae s0 a, | 
aja, = -(Qo- n) - 2 BP - n) + Z~ O- 0) 
+ p(Q, - > serge 1 (per .25) 


From equation (III.25) we can see that a AG is 
negative as long as Q, > 2n. Therefore even though 
cur calculation for aA, is crude it is negative as 
we expected by our physical argument. Since a A, 
is negative it can be seen by looking at equation 
(ETC 21) that aA, is positive and a A, is negative. 

We can rewrite the radial density equation 
({Ticto) using equactrons (Lit. 25) and (11tbsZzy) re- 


membering that our solution now is very approximate. 


Equation (III.15) becomes 


PQ PQ Q 
a(x) = Q (1 - G2 x*) - (Q- n) - G2 GE - 0) 


rQ }oy | 
+ (52 (Q.- n) - =F (Q- nD) x? (III. 26) 
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In obtaining (III.26) we have expanded the 
Bennett solution (equation (III.10)) in a binominal 
expansion and throughout the equation we kept terms 


only up tejaspowerxof x?, 


rri.4 = Calculation vorsthe: Radial) Component of the 


Internal Electric Field 


The radial component of the internal electric 
field E.. can be determined from the condition that 
in the steady state the radial flux of both type 


of charges must be equal 


ae 8 ° Cl Erec 7) 


Now, it can be seen from equation (III.6) that 


aos aie) | ae Tite2s8 
“a = + oe - HWonB, (reo, ~ De ae ( 2G) 
Therefore using equations (III.27) and (III.28) we 


obtain for the radial component of the internal elec- 
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Now rewriting equation (III.29) in terms of 


the variable x we obtain 


EL (x) = Eor (*) - Sa Ge : Gls. 3/1) 
where 
-, ia 

Eor ‘*) = Q q dx (Git... 32) 
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a is the same as in (III.30). 
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CHAPTER? LV 


DISCUSSION OF RESULTS 


Equation (III.26) can be rewritten keeping only 
the most dominant terms as 


rQ rQ 
a(x) ~» Q)(1- G2 x*) - (Q)- n) + S2(Q.- n)x*. (Ve) 


This can be rewritten as 
Q 
q (x) ~ sr Oa (> ms n) x “ GEV cz?) 


It can be seen by looking at equation (IV.2) 
that the effect of the perturbed density term is to 
reduce the density at the cylindrical axis (x = 0) 
from the Bennett value of On which is greater than 2 
to a different value n which is greater than l. 

For the region near the cylindrical axis the 
density will increase quadratically provided Oe oie 
which we have assumed. However, it should not be 
assumed that the density profile will continue to 
increase quadratically as x becomes larger. In fact 
it should not be assumed that the density profile 
will continue to increase at all for larger values 


of x. As x increases the effect of the terms which 
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contain the higher powers of x will dominate the 
solution. The effect of these terms has not been 
Kept in) Our: approximate=equacron (itr. zo0). fn Lact 

it is expected that the higher power terms will 
dominate over the quadratic terms and cause the 

radial density to decrease to a value close to the 
density at the wall given by the Bennett solution. 
Therefore, a density profile similar to that obtained 
by Daugherty and Ventrice (fig. (7)) will be obtained 
except that the density at the wall will not be zero. 
However, it is difficult to determine exactly how many 
terms in the series must be kept before this is possi- 
ble. 

Therefore, all that can be said about our appro- 
ximate solution is that it can give a radial form 
similar to the graph obtained by Daugherty and Ventrice 
(tagew(/)) for smallsyvalues oOf-rs “For-a pure m = 2 
mode to fully develop in the tube before a m = 1 mode 
a form very similar to ours is needed. To obtain a 
complete radial density profile for our problem one 
would have to include much higher powers of x in our 
solution for a A, and throughout the paper. This could 
not be done very easily analytically but could probably 


be done by computer. 
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